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Abstract 

There are several efficient methods to solve linear interval poly¬ 
nomial systems in the context of interval computations, however, the 
general case of interval polynomial systems is not yet covered as well. 
In this paper we introduce a new elimination method to solve and 
analyse interval polynomial systems, in general case. This method 
is based on computational algebraic geometry concepts such as poly¬ 
nomial ideals and Grobner basis computation. Specially, we use the 
comprehensive Gobner system concept to keep the dependencies be¬ 
tween interval coefficients. At the end of paper, we will state some 
applications of our method to evaluate its performance. 


1 Introduction 

Many computational problems arising from applied sciences deal with floating¬ 
point computation and so require to import polynomial equations containing 
error terms in computers. This redounds polynomial equations to appear 
with perturbed coefficients i.e the coefficients range in specific intervals and 
so these are called interval polynomial equations. Interval polynomial equa¬ 
tions come naturally from several problems in engineering sciences such as 
control theory [35,36], dynamical systems [30] and so on. One of the most 
important problems in the context of interval polynomial equations is to 
analyse and study the stability and solutions of an (or a system of) interval 
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polynomial(s). More generally, the problem is to carry as much as possible 
information out from an interval polynomial system. Many scientific works 
are done in this direction using interval arithmetic [2], for instance compu¬ 
tation of the roots in certain cases [6,7], however they do not enable us to 
obtain the desired roots, at least approximately [6]. Another example consists 
those works which contain (the most popular) method to solve an interval 
polynomial equation by computing the roots of some exact algebraic poly¬ 
nomials, while it is hard to solve an algebraic equation of high degree which 
has its own complexity challenging problems [9,11,12,18,21], There is also a 
new method described in [41] which counts the zeros of a univariate inter¬ 
val polynomial. In addition to numerical methods, there are some attempts 
to combine numeric and symbolic methods to solve an interval polynomial 
system. In [8], Falai et al. state a modification of Wu’s characteristic set 
method for interval polynomial systems, and use numerical approximation 
to find an interval containing the roots. The essential trick in this work is to 
omit all the terms with interval coefficients containing zero, which permits 
the division of interval coefficients simply. This consideration may fail some 
important polynomials. 

In this paper, we try to use exact symbolic methods to facilitate analysing 
interval polynomial systems thanks to algebraic elimination methods like 
parametric computation techniques to analyse a parametric polynomial sys¬ 
tem which allows us to consider all exact polynomials arising from an interval 
polynomial. As we will state later, it is very important to keep the trace of 
interval coefficients during the computations. Roughly speaking we asso¬ 
ciate an auxiliary parameter to each interval coefficient provided that each 
parameter ranges over its own related interval only. Nowadays there are 
important results [38,39], efficient algorithms [17,19,20,22,24,25] and power¬ 
ful implementations [33,34] in the context of parametric computations and 
analysing parametric polynomial systems. We introduce the new concept 
interval Grobner system, for a system of interval polynomials using the con¬ 
cept of comprehensive Grobner systems [39] which is used to describe all 
different behaviours of a parametric polynomial system. Interval Grobner 
system contains a finite number of systems where each one is a Grobner 
basis for (non-interval) polynomial systems obtained from the main system. 
It is worth noting that against to [8], we don’t omit any interval coefficient 
and cover all possible cases for the exact coefficients arising from the inter¬ 
vals. We also design an algorithm to compute interval Grobner system of an 
interval polynomial system. 
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This paper is organized as follows. In Section 2 we state introductory 
definitions and recall interval arithmetic. In Section 3 we start to explain 
interval polynomials and their related concepts. We next receive to Section 4 
which states the main idea behind this paper. To recall the concepts of com¬ 
putational algebraic tools we state Section 5 containing a brief introduction 
to Grdbner basis and comprehensive Grobner system, together with their 
related algorithms. After, we describe our elimination method for interval 
polynomial systems in Section 6. Finally, as some applications and example, 
we state the Section 7 which contains two applied examples. 


2 Preliminaries 


In this section we recall the interval arithmetic and related concepts which 
are needed for the rest of this text. The main references of this section are 
[16] and [26]. Let M denote the set of real numbers while M* is used to show 
the extended real numbers set i.e. M U {—oo, oo}. 


Definition 2.1. Let a,b E M*. We define 4 kinds of real intervals defined by 
a and b as follows: 


Closed interval 
Left half open interval 
Right half open interval 
Open interval 


[a, b] — {x | a < x < b} (a, b ±oo) 

(a,b] = {x \ a < x < b} (b ±oo) 

[a, b) = {x \ a < x < b} (a / ±oo) 

(a, b) = {x | a < x < b} 


It is worth noting that approximately all of existing texts in the subject 
of interval computation deal with closed intervals. Most of times we denote 
the intervals by capitals, and their lower (resp. upper) bounds by underbar 
(resp. overbar), as 

X = [X,X] 

to denote closed intervals. However, as there are some applied problems 
including non-closed intervals we preferred to consider all different types of 
intervals. 


Remark 2.2. Having all different kinds of intervals at once, we use the 
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notion [a,b,i,j] where i,j G {0,1} to denote the intervals in (1), as follows: 


a,b,i,j} 


[a, b] if i — j — 1 

(a, b] if i — 0,j — 1 

[a, b) if i — 1 , j — 0 

(a, 6 ) if i = j = 0 


However when all of intervals come from one sort of presentation, we prefer 
to use the (1) form. 

Now we recall the interval arithmetic and discuss on the interval depen¬ 
dencies what will occur in the evaluation of interval expressions. 

There exist two equivalent ways to state interval arithmetic. The first is 
based on the endpoints of intervals while the second considers each interval 
as a subset of real numbers. Let [ai, b\, and [ 02 , b 2 , * 2 ,^ 2 ] be two real 
intervals. Note that each real number a is considered as [a, a, 1,1] which is 
called a degenerate interval. Four essential arithmetic operations are defined 
as follows: 


[°i) bi,ii,ji] + [a 2 ,b 2 ,i 2 , j 2 \ 
[ai,bi,ii,ji\ — [a 2 , b 2 , 12 , J 2 ] 
[ai,bi,ii,ji] x [a 2 ,b 2 ,i 2 , J 2 ] 


[ai + a 2 , h + b 2 , min(i 1 , i 2 ), min (j 1 ,j 2 )] 
[ai - b 2 , bi - a 2 , min(ii, j 2 ), min(ji, i 2 )\ 
[a k b e , a k >be>, min (i k , j e ), min (i k >,je')} 


where a k bi and a k /b^ are the minimum and maximum of the set {aia2, a\b 2 , b\a 2 , 61^2} 
respectively, and finally 

[a 1 ,b 1 ,i 1 ,j 1 \/[a 2 ,b 2 ,i 2 ,j 2 \ = [oq, b x , h, ji] x [l/6 2 , l/a 2 , J2, *2] 


provided that a 2 > 0 or b 2 < 0 or a 2 = i 2 = 0 or b 2 = j 2 = 0 . 

Note in the above relations that all ambiguous cases 00 — 00 , ±00 x 0, 
and jj will induce the biggest possible interval i.e. M. 

As an easy observation, when an interval X = [a,b,i,j] with a 7 ^ b 
contains zero, we can compute 1/A" as follows: 


• If a = 0 then 


1 1 
A" [a, b, 0, j] 


[5- 


+ 00 , j, 0], 
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• If b = 0 then 


1 


1 

X 


r , • n1 = [-00, —, 0, ®], 

a, b, i, (JJ a 


• If ab < 0 then by seperating X as X = [a, 0, i, 0] U [0, b, 0, j] we have 


1 1 

X [a, 0, i, 0] U [0, b, 0, j] 


r 1 , .1 

-oo, —, 0, *J U 7 ,+oo, j,0 . 
a b 


Now consider A and B are two intervals as two sets of real numbers. We can 
state the above definitions of four essential arithmetic operations as 

A op B = {x | 3 a G A, b e B : x — a op b} 

where op € {+, —, x, /}. 

Remark 2.3. Although interval arithmetic seems to be compatible with real 
numbers arithmetic, but this affects distributivity of multiplication over addi¬ 
tion and the existance of inverse elements. More preciesly for each intervals 
X, Y and Z, 

• X x (Y + Z) C X x Y + X x Z 7 
and if X is non-degenerated then 


• I x I / 1, but 1 G X x y and 


• X + (-X) ^ 0 ; but 0 G X + {-X). 


Furthermore, if X contains some negative real numbers, then 

X n ^X_^X. 

n times 

To see this, let for instance X = [a,b,i,j\ where a < 0 and |a| < b. Then, 
X 2 = [0, b 2 ,l,j\ while X x X = [ab, b 2 , min(i, j), j]. To solve this inconsis¬ 
tency, we define the wth power of an interval for non-negative integer n as 
follows: 


a,b,i,j] 


n 


1 n — 0 

[a n ,b n ,i,j] 0 <a 

[b n ,a n ,j, i] b< 0 

[ 0 , max(a", b n ),l,c\ a < 0 < b 


where c 


i | 6 | < |a| 
j otherwise 
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Let us now evaluate some expressions to illustrate more challenging prob¬ 
lems dealing with interval arithmetic. Let f(x,y) = X = [1,2] and 
Y = [1,3]. We compute f(X,Y) in two ways. The first way is to compute 
/( X,Y) as an usual evaluation using interval arithmetic: 

jL y = [1.2]/[2,5] = [1/5,1], (2) 

However, one can manipulate the expression to see 

—Lp = Yyx = V P/ 2 - 4 ) = [!/4,2/3] (3) 

Let us separate x and y first depending on f(x,y ): we call y (resp. x), a 
first (resp. second ) class variable of / as it appears one (resp. more than 
one) time in the structure of /. As it can be easily seen, the answer of (3) 
is a narrower interval and in fact the exact value. The reason is that X is a 
second class variable for (2) and so it brings dependency between two parts 
of the expression. This is while there is no dependency in (3) given that A 
appears only one time, and so it is a first class variable. Dependency is one of 
the crucial points of this paper to find the solution set of interval polynomial 
systems. 

Dependencies are the main reason that causes to appear an amount of 
error by introducing larger intervals than the exact solution as all publications 
in this area try to avoid dependencies. Nevertheless it is possible to cancel 
dependencies by considering A" — A" = 0 (see [16]) as well as X/X = 1 easily, 
while sometimes this becomes a difficult work. 


3 Interval Polynomials 

Let R be the field of real numbers, considered as the ground field of compu¬ 
tations all over the current text and consider the set {xj,... ,x n } to be the 
set of variables. 

Definition 3.1. Each polynomial of the form 

m 

[/] = 53 °*’ hi ' ' ' ' X n n ( 4 ) 

2=1 
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is called an interval polynomial, where [a*, hi, kj\ is a real interval for each 
i — 1 ,..., m, and each power product xfi 1 ■ ■ ■ xlf in is called a monomial where 
the powers are non negative integers. We denote the set of all interval poly¬ 
nomials by [R] [xi,..., x n ]. 

Definition 3.2. Let [/] he an interval polynomial as defined in (f). The set of 
all polynomials arising from [/] for different values of intervals in coefficients 
is called the family of [/] and is denoted by tF([f]). More preciesly: 

m 

F{[f]) = {V, Cjxf 11 ■ ■ ■ x°fi n I a G [ai^^kifi = 1 ,..., m} 

i— 1 

Similar to the family of an interval polynomial, we can define the family 
of a set of interval polynomials as follows: 

Definition 3.3. Let S := {[/]i, ■ • ■, [/]«} he a set of interval polynomials 
with fF{[f]j) = Tj for each j — 1,..., i. We define the family of S to be the 
set T\ x • • • x Tt, denoted by IF(S). 

We now define the concept of solution set for an interval polynomial. 

Definition 3.4. For an interval polynomial [f] G [R] [x \,..., x n \, we say that 
a = (ai,..., a n ) G R n is a real solution or a real root of [/], if there exists a 
polynomial p G -R([/]) such that p(a) = 0. Similarily, we say that a system 
of interval polynomials has a solution, if the contained interval polynomials 
have a common solution. 

Example 3.5. Let us find the solution set of 

[—2, —l]x 2 + [1, 2]x + [1, 3] = 0 

where all intervals are closed. When x > 0, we have 

[-2, -l]x 2 + [1, 2]x + [1, 3] = [—2x 2 + x + l,-x 2 + 2x + 3] = [0, 0] 

So we must have 


0 < — 2x 2 + x + 1 and — x 2 + 2x + 3 < 0 


which implies 


x G [1, 3]. 
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Similarly when x < 0, we have 

[-2, -l]x 2 + [1, 2]x + [1, 3] = [ -2x 2 + 2x + l,-x 2 + x + 3] = [0, 0] 
or equivalently 


0 < —2a ; 2 + 2a; + 1 and — a ; 2 + a; + 3<0 


which concludes 

x G [-1.31,-0.36], 

Thus the solution set of this interval polynomial is 


-1.31,-0.36] U [1,3]. 


It is notable that using interval arithmetic in the well-known solution way of 
a quadratic polynomial equation due to discriminant, we receive to 

[-2.15,-0.06] U [0.382,8.4] 

as the solution set that contains an amount of error. 

4 The idea 

In this section we are going to describe the problems which may occur using 
the usual elimination method on a system of interval polynomials. To facili¬ 
tate the description, let us give an example. Consider the system containing 

fi = [1, 2]x\ + x 2 + 2 x 3 , /2 = [1,4]a;i + x 2 + l, fa = [3,4]xi + x 2 + 4a; 3 . 

Going to eliminate the variable x 2 , we have 

f\ ~ f '2 = [-3, l]xi + 2a ; 3 - 1, / 3 — / 2 = [-1, 3]ici + 4a ; 3 - 1. 

One may now conclude that if we choose 0 from both intervals [—3,1] and 
[—1,3] then the system has no solution because 2a ; 3 — 1 = 4a ; 3 — 1 = 0. 
However this case is impossible since both of intervals [—3,1] and [—1, 3] can 
not be zero at once! To see this, notice that [—3,1] comes from [1,2] — [1,4] 
and so for [—3,1] being zero, [1,4] must give some values in [1,2]. On the 
other hand, [—1,3] comes from [3,4] — [1,4] and so this interval can be zero 
only when [1,4] gives some values in [3,4] and this is a contradiction. This 



simple linear system shows that the usual elimination method can cause a 
wrong decision or appearing some extra values in the solution set. The main 
reason is that we forgot the dependencies between [—3,1] and [—1,3] during 
the computation while they are both dependent on [1,4] and so they are 
dependent. 

To solve this problem, we must keep the trace of each interval coefficient. 
In doing so, our idea is to use a parameter instead of each interval, to see 
how new coefficients are built. For instance let us substitute [1,2], [1,4] and 
[3,4] by a, b and c as parameters in the above example. So we have 

/i = ax i +x 2 + 2 x 3 , f 2 = bx i +x 2 + l, f 3 = cx i + x 2 + 4x 3 
and doing elimination steps we have 

h - h = (a - b)x i + 2 x 3 - 1 , f 3 - f 2 = (c - b)x i + Ax 3 - 1 . 

Now one can conclude that under the assumption that a — b = 0, the coeffi¬ 
cient c — b can not be zero. The reason is that if c — b = 0 then a = c while 
1 < a < 2 and 3 < c < 4. Therefore using parameters prevent us taking 
wrong decisions about the solution set. 

The main question here is that how we can use elimination method when 
the coefficients contain some parameters. In fact as we will state in the next 
sections, we do the elimination steps thanks to Grobner basis and for the 
parametric case, we use the concept of comprehensive Grobner system, to 
see the simplest possible polynomials to solve. So our idea is to convert 
the interval polynomial systems to a parametric polynomial system and use 
the parametric algorithmic aspects, of course with some modifications, to 
solve the parametric system by dividing the solution set into finitely many 
components. At the end, we convert the result to see the solution set of the 
interval polynomial system. 


5 Grobner Bases and Comprehensive Grobner 
Systems 

In this section we recall the concepts and notations of ordinary and para¬ 
metric polynomial rings. Let K be a held and x \,..., x n be n (algebraically 
independent) variables. Each power product x “ 1 • • • x^ n is called a monomial 
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where aq ,a n E Z> 0 . Because of simplicity, we abbreviate such monomi¬ 
als by x“ where x is used for the sequence xi,... ,x n and a = (aq,..., a n ). 
We can sort the set of all monomials over K by special types of total orderings 
so called monomial orderings, recalled in the following definition. 

Definition 5.1. The total ordering -< on the set of monomials is called a 
monomial ordering whenever for each monomials x",x i3 and x 7 we have: 

• x“ -< x^ x 7 x" -< x 7 x /3 , and 

• -< is well-ordering. 

There are infinitely many monomial orderings, each one is convenient for 
a special type of problems. Among them, we point to pure and graded reverse 
lexicographic orderings denoted by -<i ex and -< gre viex as follows, assume that 
x n -<■■■-< x\. We say that 

• x a -<i ex x 5 whenever 


ai = /?i, ■ ■ ■, oti = /% and a i+ i < fi l+ \ 


for an integer 1 < i < n. 


~^grevlex 


iP if 


n 


n 


y a% < 

i= 1 i —1 

breaking ties when there exists an integer 1 < i < n such that 


Pri •> • • • i &n—i fin—i arid Oi n _{_i j. 


It is worth noting that the former has many theoretical importance while the 
latter speeds up the computations and carries fewer information out. 

Each K— linear combination of monomials is called a polynomial on aq,..., x n 
over K. The set of all polynomials has the ring structure with usual poly¬ 
nomial addition and multiplication, and is called the polynomial ring on 
xi,... ,x n over K and denoted by K[x i, ..., x n ] or just by A'[x], Let / be 
a polynomial and -< be a monomial ordering. The greatest monomial w.r.t. 

-< contained in / is called the leading monomial of /, denoted by LM(/) 
and the coefficient of LM(/) is called the leading coefficient of / which is 
pointed by LC(/). Further, if A is a set of polynomials, LM(A) is defined 
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to be {LM(/)|/ G F} and if / is an ideal, in(I) is the ideal generated by 
LM(/) and is called the initial ideal of I. We are now going to remind the 
concept of Grobner basis of a polynomial ideal which carries lots of useful 
information out about the ideal. 

Definition 5.2. Let I be a polynomial ideal of K [x] and -< be a monomial 
ordering. The finite set G C I is called a Grobner basis of I if for each non 
zero polynomial f G I, LM(/) is divisible by LM(g) for some g E G. 

Using the well-known Hilbert basis theorem (See [4] for example), it is 
proved that each polynomial ideal possesses a Grobner basis with respect 
to each monomial ordering. There are efficient algorithms also to compute 
Grobner basis. The first and the most simplest one is the Buchberger al¬ 
gorithm which is devoted in the same time of introduction of Grobner basis 
concept while he most efficient known algorithm is the Faugere’s F 5 algorithm 
[10] and another signature-based algorithms such as G 2 V [13] and GVW [14], 
It is worth noting that Grbbner basis of an ideal is not unique necessarily. To 
have unicity, we define the reduced Grobner basis concept. As an important 
fact the reduced Grobner basis of an ideal is unique up to the monomial 
ordering. 

Definition 5.3. Let G be a Grobner basis for the ideal I w.r.t. -<. Then 
G is called a reduced Grobner basis of I whenever each g E G is monic, i.e. 
LC(g) = 1 and none of the monomials appearing in g is divisible by LM(h) 
for each h G G\ {g}. 

One of the most important applications of Grobner basis is its help to 
solve a polynomial system. Let 

[ h = o 
{ fk = o 

be a polynomial system and / = (fi,---,fk) be the ideal generated by 
/i,..., fk- We define the affine variety associated to the above system or 
equivalently to the ideal / to be 

V(/) = V(A, ...>/*) = {a G in/i(aO = • • • = /*(a) = 0} 

where K is used to denote the algebraic closure of K. Now let G be a 
Grobner basis for / with respect to an arbitrary monomial ordering. As an 
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interesting fact, / = ( G ) which implies that V(/) = V(G). This is the key 
computational trick to solve a polynomial system. Let us continue by an 
example. 

Example 5.4. We are going to solve the following polynomial system: 

{ ,x 2 — xyz + 1 = 0 
y 3 + z 2 - 1 = 

xy 2 + z 2 =0 

By the nice properties of pure lexicographical ordering, the reduced Grobner 
basis of the ideal I = (x 2 — xyz + 1, y 3 + z 2 — 1, xy 2 + z 2 ) C Q[x, y, z] has the 
form 

G = { 9 i(z),x ~ g 2 (z),y - g 3 (z)} 
w.r.t. z -< iex y -< tex x, where 

gi(z) = z 15 — 3 z 14 + 5 z 12 — 3z 10 — z 9 — z 8 + 4 z 6 — 6 z 4 + 4 z 2 — 1 
g 2 {z) = 2z 14 - 9z 13 + llz 12 + 2z n - 7z 10 - 3z 9 + 2z 8 - z 7 + 4z 6 + 

+ 7 z 5 - 10z 4 - 6z 3 + llz 2 + 2z - 4 

_ g 3 (z) = z 13 -3z 12 + z 11 + 2z 10 + z 9 -z 8 -2z 6 + 2z 4 -z 3 -3z 2 + 1. 

This special form of Grobner basis for this system allows us to find V( G) by 
solving only one univariate polynomial gi(z) and putting the roots into the 
two last polynomials in G. 

Suppose now that the same system of Example 5.4 is given as follows 
with parametric coefficients on parameters are a, b and c: 

{ ax 2 — ( a 2 — b + 1 )xyz + 1 = 0 
y 3 + c 2 z 2 - 1 = 

(a + b + c)xy 2 + z 2 =0 

The solutions of this system depend on the values of parameters apparently as 
we can see that the system has no solutions whenever a = 0 and 6=1 while 
it converts to the system of Example 5.4 for a = 1,6 = 1 and c = — 1 and 
so has some solutions. To manage all of different behaviours of parameters 
which cause to different behaviour of the main system, we recall the concept 
of comprehensive Grobner system in the sequel. By this we can divide the 
space of parameters, i.e. K into a finite number of partitions, for which the 
general form of polynomials contained in assigned Grobner basis is known. 
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Let K be a field and a := cq,..., a t and x := xi ,..., x n be the sequences 
of parameters and variables respectively. We call K [a] [x], the parametric 
polynomial ring over K, with parameters a and variables x. This ring is in 
fact the set of all parametric polynomials as 

m 

i= 1 

where p, G A'[a] is a polynomial on a with coefficients in K, for each i. 

Definition 5.5. Let I C K [a] [x] be a parametric ideal and -< be a monomial 
ordering on x. Then the set 

G(I) = {(Ei, Ni,Gi)\i = l,...,e}<z K[ a] x K[ a] x K[ a][x] 

is said a comprehensive Grobner system for I if for each (Ai,..., A*) G K* 
and each specialization 

^(A..a*) : Af[a][x] -> W[x] 

m 

E!:=i^ Pi ( Xl ’' ' ' ’ 

i= 1 

there exists an 1 < i < l such that (Ai, ...,A t ) G ’V(E i ) \ V(W) and 
ct(Ai,...,a t )(Gj) is a Grobner basis for c r (Ai,...,A t )(-^) with respect to -<. Because 
of simplicity, we call Ei and Ni the null and non-null conditions respectively. 

Remark that, by [39, Theorem 2.7], every parametric ideal has a compre¬ 
hensive Grobner system. Now we give an example from [20] to illustrate the 
definition of comprehensive Grobner system. 

Example 5.6. Consider the following parametric polynomial system in Q[a, b, c] [x , y]: 

{ ax — b =0 
by -a =0 
cx 2 — y = 0 
cy 2 — x = 0 

Choosing the graded reverse lexicographical ordering y -< x, we have the fol¬ 
lowing comprehensive Grobner system: 
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Gi 

Ei 

Ni 

{1} 

{} 


{bx — acy, by — a} 


m 

{cx z — y, cy z — x} 

{a,b} 

w 

{x,y} 

{a,b, c} 

ii 


For instance, for the specialization (7(i,i,i) for which a H > 1, b (->• 1 and 
c i —y 1, 

ffii,i,i)({&c - ac V, - a}) = {x - y, y - 1} 
is a Grobner basis of <7(1,1,1) ((£)). 

It is worth noting that if \{Ef) \ V(IVj) = 0 for some i, then the triple 
(. Ei , iVj, Gj) is useless, and so it must be omitted from the comprehensive 
Grobner system. In this case we say that the pair (E^, Nf) is inconsistent. It is 
easy to see that inconsistency occurs if and only if N { c \J(Ef) and so we need 
to an efficient radical membership test to determine inconsistencies. In [19, 
20] there is a new and efficient algorithm to compute comprehensive Grobner 
system of a parametric polynomial ideal which uses a new and powerful 
radical membership criterion. Therefore we prefer to employ this algorithm 
so called PGB algorithm in our computations. Another essential trick which 
is used in [20] is the usage of minimal Dickson basis which reduces the content 
of computations in PGB. Before explain it, let us recall some notations which 
are used in the structure of PGB. Let -< x and -< a be two monomial orderings 
on K [x] and K [a] respectively. Let also -< x , a be the block ordering of -< x 
and -< a , comparing two parametric monomials by -< x , breaking tie by -< a . 
For a parametric polynomial / G K [a] [x] , we denote by LM X (/) (resp. by 
LC X (/)) the leading monomial (resp. the leading coefficient) of / when it is 
considered as a polynomial in Jl(a)[x], and so LC X (/) G K[ a]. 

Definition 5.7. By the above notations, let P C K [a] [x] be a set of para¬ 
metric polynomials and G C P. Then, G is called a minimal Dickson basis 
of P denoted by MDBasis(P), if: 

• For each p G P, there exist some g G G such that LM x (g) | LM x (p) 
and 

• For each two distinct polynomials in G as g\ and g 2 , none o/LM x ((? 1 ) 
and LM x (g 2 ) divides another. 
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The case which occurs in PGB to compute a minimal Dickson basis for P 
is only when P is a Grobner basis for ( P ) itself w.r.t. -< x a and PnK[a] = {0}. 
In this situation, it suffices by Definition 5.7 to omitt all polynomials p from 
P for which there exists a p' E P such that LM x (p') | LM x (p). 

The PGB algorithm as is shown below, uses PGB-MAIN algorithm to 
introduce new branches in computations. 


Algorithm 1 PGB 
1: procedure PGB(P, ~<a, ~<x) 

2: E,N:={ 

3: -<x,a : =The block ordering of -< x , -< a 

4: Return PGB-main(P, E, N, -< x ,a); 

5: end procedure 


The main work of PGB-main is to create all necessary branches and 
import them in comprehensive Grobner system at output. In this algorithm 
A* B is defined to be the set {ab \ a e A, b e B}. 
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Algorithm 2 PGB-main 
procedure PGB-MAIN(P, E, N, -< x ,a) 

G := The reduced Grobner basis for P w.r.t. -< ax 
if 1 G G then 

Return (E,N, {1}); 
end if 

G r :=G n K[ a]; 

if IsConsistent(P, N * G r ) then 
PGB:={(E,N*G r ,{ 1})}; 
else 

PGB := 0; 

end if 

if IsCoNSlSTENT(G r , N ) then 
G m :=MDBasis(G \ G r ); 
else 

Return (PGB); 
end if 

h := lcm(h],..., hk), where h t = LC x (p,) and g t G G' m ; 
if IsCoNSlSTENT(G r , N * {/i}) then 
PGB := PGB U {(G r , N * {h}, G m )}; 

end if 

for i = 1,..., k do 

PGB := PGB U PGB-main(G\G>, G r U{hi}, N* {Jl'I 1 , hj], ^ a ,x) 

end for 

end procedure 


As it is shown in the algorithm, it computes first a Grobner basis of 
the ideal ( P) over K[a., x] i.e. G, before performing any branches based on 
parametric constraints, according to [20, Lemma 32] as follows: 

Lemma 5.8. By the notations used in the algorithm, for each specialization 
(A,,...,A,)EV(GJ\V( JJ LC x (</)) 

g£G\G r 

then cr(Aj ,...,\ t )(G) is a Grobner basis for a^\ lt „^\ t )({P)). 

After this, the algorithm computes a minimal Dickson basis i.e. G m and 
continues by taking a decision for each situation that one of the leading 
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coeffiecients of G m is zero. By this, PGB-main constructs all necessary 
branches to import in comprehensive Grobner system. All over the algo¬ 
rithm, when it needs to add a new branch (TJj, iVj,Gj) into the system, the 
algorithm IsCONSISTENT is used as follow to test the consistency of para¬ 
metric conditions ( E{,Ni ). 


Algorithm 3 IsConsistent 
procedure IsConsistent(17, N ) 
flag :=false; 

for g E N while flag =false do 
if g f a/ (E) then 
flag :=true; 

end if 
end for 
Return flag ; 
end procedure 


The main part of this algorithm is radical membership test. The powerful 
trick which is used in [19,20] to radical membership check is based on linear 
algebra methods tackling with a probabilistic check. We refer the reader to 
[20, Section 5] for more details. 


6 Elimination Method to Solve Interval Poly¬ 
nomial Systems 

In this section we introduce the new concept of interval Grobner system and 
its related definitions and statements. 

Now we state the following proposition as an immediate consequence of 
Definition 3.4. Recall that for a polynomial system S C Mfoq,..., x n ] the 
variety of S is the set of all complex solutions of S, denoted by V(S). 

Proposition 6.1. A system [S'] of interval polynomials has a solution if and 
only if there exists a polynomial system S in ,F([S]) with V(S) 0. 

There is an efficient criterion due to the well-known Hilbert Nullestelen- 
satz theorem which determines if V(S) 0 by Grobner basis: V(S) 0 if 

and only if the Grobner basis of (S) does not contain any constant. Note 
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that there are infinitely many polynomial systems in X([S]) for an interval 
polynomial system [S] and so it is practical impossible to check all of them by 
Nullestelensatz theorem. Nevertheless, we give a finite partition on the set of 
all polynomial systems arising from [S] using the concept of comprehensive 
Grobner system. 

Definition 6.2. Let [S'] = {[/]i,..., [f]e} be a system of interval polynomi¬ 
als. We define the ideal family of [.S'], denoted byllF([S}) to be the set 

IT(\S]) = {(p 1 ,...,p l ) | (p u ... ,pi) e ^([S])} 

Theorem 6.3. Let [S] be a system of interval polynomials and -< be a mono¬ 
mial ordering on M[xi,..., x n ]. Then 

• The set of initial ideals { in(I ) | / e XX([S])} is a finite set, and 

• For each set of ideals oflE([S}) with the same initial ideal, there exists 
a set of parametric polynomials which induces the ideals by different 
specializations. 

Proof. To prove this theorem, we use the concept of comprehensive Grobner 
system. Suppose that S* is obtained by replacing each interval coefficient 
by a parameter. Note that if an interval appears in t > 1 coefficients, then 
we assign t distinct parameters to it. It is easy to check that each element 
of XX([S]) is the image of S* under a suitable specialization. On the other 
hand by [39, Theorem 2.7] S* has a finite comprehensive Grobner system as 
Q = {(Ei , Ni, G i),..., (E k , N k , Gk)}, where for each specialization a there 
exists a 1 < j < k such that LM(cr(S*)) = LM(Gj). It is worth noting that 
although there is a finite number of branches in Q, we can also remove those 
specializations with complex values, and also those with values out of the 
assigned interval. Thus for each / e XX"([S']) there exists an 1 < i < k with 
in(I ) = (LM(Gj)) and this finishes the proof. □ 

What is explained in the proof of Theorem 6.3 yields to extend the concept 
of comprehensive Grobner system for interval polynomials. 

Definition 6.4. Let [S] C [M] [x i,..., x n ] be a system of interval polynomials 
with t interval coefficients, and -< be a monomial ordering on Mfxi,..., x n ]. 
Let also that Q = {( Ei , Ni,G i),..., (E k , N k , G k )} be a set of triples 

(Ei, N t , Gf) G R[h u ...,h t \ x R[h u ...,h t \ x R[h u ..., h t \[x i,... ,x n ] 
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where {h\,... ,h t } is the set of parameters assigned to each interval coeffi¬ 
cient. Then we call Q an interval Grobner system for [S'] denoted by ^([S 1 ]) 
if for each t—tuple (ai,..., af) of the inner values of interval coefficients there 
exists an 1 < i < k such that: 


• For each p G E % , p{a \,..., a t ) = 0, 

• There exist some q G iV* such that q(a\, ■ ■ ■ ,a t ) ^ 0, and 

• cr(Gi) is a Grobner basis for (cr([S'])) with respect to -<, where a is the 
specialization hj (->• aj for j = 1 ,... ,t. 

Theorem 6.5. Each interval polynomial system possesses an interval Grobner 
system. 

Proof. Let S* be the parametric polynomial system obtained by assigning 
each interval coefficient to a parameter. As mentioned in the proof of The¬ 
orem 6.3, ^([5]) is the same comprehensive Grobner system of S* where 
each parameter is bounded to give values from its assigned ideal. On the 
other hand it is proved that each system of parametric polynomials has a 
comprehensive Grobner system, which terminates the proof. □ 


We give now an easy example to illustrate what described above. 


Example 6.6. Consider the interval polynomial system 


[ 5 ] 


[—1, 2)xy + [0, l)y + [3,5) = 0 
[—3, l)xy 2 + [1, 3)y = 0 


(5) 


To obtain a parametric polynomial system, we assign the intervals [—1, 2), [0,1), [3, 5), [—3,1) 
and [1, 3) by hi,..., h 5 respectively. Then we observe the parametric polyno¬ 
mial system 


S* = {hixy + h 2 y + h 3 , h A xy 2 + h 5 y} Cl[/ii,..., h 5 ][x, y] 


Using the lexicographic monomial ordering y -< x we can compute a compre¬ 
hensive Grobner system for ( S*) which contains about 19 triples. However 
some of them are admissible only for some values of parameters out of their 
assigned interval. For instance the triple ({1}, {hi, h 2 , h^, h$}, {h 3 }) is not 
acceptable in this example, since h§ G [1,3] and so it can not be zero. By 
removing such triples, there remains only 8 one shown in the following table. 
Therefore the following table shows £L;([S']). 
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Ei 

Ni 

Gi 

{h 3 h 5 } 

{h\, h 2 , h 4 } 

Ul 

{h 3 h 4 } 

{h \, h 2 } 

{1} 

{hih 2 h 3 h 4 } 

{h 3 h 4 — h 4 h§} 

Ul 

{h2h 3 h 3 } 

{hi,h 4 } 

{1} 

{h\h 3 h 3 } 

{h 4 } 

{1} 

{h 2 h 3 h 4 } 

{hi} 

{h 2 y + h 3 , h^h 4 x - h 2 h 5 h 3 } 

{h\h 4 } 

{ h - 2 , h 3 h 4 — h\h 3 } 

{h 3 + hixy} 

{hih 2 h 3 h 4 {h 3 h 4 — h\h 3 )} 

{} 

{h\xy + h 2 y + h 3 ,-h 4 h 2 y z - h 4 h 3 y + h 5 h 4 y, 
—h 4 h\h 4 x + h\h§h 3 x + h 2 h 4 h 3 y + h 2 h 4 lif i , 
—h 4 h 2 h 3 y - h\h 4 + h 4 h 5 h 3 } 


Table 1: Interval Grobner system of System (5) 

6.1 Computing Interval Grobner Systems 

In tis section we state our algorithm so called IGS to compute interval 
Grobner system for an interval polynomial system. This algorithm is based 
on the PGB algorithm with some extra conditions for the definition of con¬ 
sistency. To begin let [S] = {[/]i,..., [f]e} C [R] [x \,..., x n ] be a system of 
interval polynomials, where for each 1 < j < £ } 

mj 

[fh = 

i =1 

and (afij-,1, • • •, € Z" 0 , for each i. As it is mentioned in Theorem 6.3, we 

assign to each interval coefficient [a. VJ , b t] \ a parameter h t j to convert [S'] to 
a parametric polynomial system S*. The following proposition describes the 
relations between comprehensive Grobner systems of S* and interval Grobner 
bases of [S]. 

Proposition 6.7. Using the above notations, let [Q] and Q be an interval 
Grobner basis for [S] and a comprehensive Grobner basis for S* respectively 
w.r.t. the same monomial ordering. Then for each ( E,N,G ) G [Q], there 
exists (£', N', G') e Q such that V(E) \ V(N) C Y(E') \ V(N') and G, G’ 
have the same initial ideal. 

Proof. This comes from Definitions 6.4 and 5.5. □ 
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According to the above proposition, to compute an interval Grobner basis 
for [S], it is enough to compute a comprehensive Grobner basis for S*, and 
use a criterion to omit those triples ( E,N,G ) lying in Q \ [Q \, we call them 
redundant triples. 

Remark 6.8. Note that for a triple ( E,N,G) in Q \ [Q], the intersection of 
V (E) \ V(. N ) with the cartesian product of interval coefficients is empty. 

We are now going to present a criterion to determine the elements of 
Q\\Q]. This criterion of course is based on the answer of this question: 

How can we sure that a system of polynomials E C R[cq,. .., a t \ has a real 
root in the interval [aq, fii) x • • • x [a t , /3 t ) ? 


In the case for which (E) is zero dimensional, this question is answered 
totally thanks to some efficient computational tools like Sturm’s chain by 
isolating the real roots. However in the case of positive dimensional, there 
exist only some algorithm to isolate the real roots. Among them there exists 
an algorithm which determines whether a multivariate polynomial system 
has real root or not. 

Because of the reasons declared above, we convert the above key question 
to the problem of determining whether a polynomial system has a real root 
or not. 


Theorem 6.9. Let E C • . • j be a finite polynomial set . Let also 
F = EU{at + (at- Pi)bf - ai \ i = 1,..., t} C R[aq,- a,.h\ r > ; .., ly] 


where bj’s are algebraic independent by at’s and [cq, (df) be a real interval for 
each i — 1,..., t. Then the system E = 0 has a solution in [aq, fdi) x • • • x 
[a*, {3 t ) if and only if the system. F = 0 has a real solution. 


Proof. Let E 
also 


0 has a solution (71 ,... ,7 1 ) € [aq,/3i) x • • • x [a t ,Pt)- Let 


Vi = 


J 


otj - 7 i 

Ti- Pi 


for each i = 1,..., t. It is easy to see that 


7 i + (7 i ~ Pi)Vi ~oti = 0 
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which implies that (71,..., 7 1 , rji ,..., rj t ) is a solution of F — 0 . 

Conversely, suppose that there exists (71,.. , y t , r/i,..., ry) G R 2t which 
is a solution of F — 0 , i. e. /(71,..., 7 1) = 0 for each / G F and 7 \ + (7* — 
Pi)rp — = 0 , for each i — 1 ..., t. It is enough to show that 7 \ G [a*, /?;). 
In doing so, we see that 


7. = 


otj + foyf 

1 + r /, 2 


(A 



+ 


2 

Indeed, 0 < < 1 and this shows that a* < (/?* — a*) 

' 'i 


Vi 


1 + vf 

% 


+ ai < /3i 


which finishes the proof. 


□ 


Remark 6.10. Note that for the intervals [a, 00) and (—oo,/3] we can use 
the auxiliary polynomials a — a — b 2 and a — /3 + b 2 respectively. 

Using Theorem 6.9 and the Remarks 6.8, 6.10 we can determine the 
elements of Q \ [Q\ exactly (see the notations of Proposition 6.7). 

Corollary 6.11. Let ( E,N,G ) G Q and [07 ,/?i),..., [at, fit) be t real inter¬ 
vals. Then (. E , N , G ) is redundant if and only if the system F = 0 has no 
real roots, where 

F = EU{ai + (at- fdfb- - <37 | i = 1,..., t} U {]^[(c fl r/ - 1)} 

g&N 

C R[a 1? ..., at, h ,..., b t , c g : g G N]. 

Proof. The proof comes from Theorem 6.9 and this fact that if Y\ g&N {c g g — 
1) = 0 then there exists a g G N for which c g g — 1 = 0 which implies that 
< 7 ^ 0 . □ 

The above corollary is the criterion which determines all redundant triples, 
and so tackling this criterion with PGB algorithm we can design our new 
algorithm to compute interval Grdbner systems. We design now the IGS 
algorithm by its main procedure. 
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Algorithm 4 IGS 
procedure IGS(S'[ ], -< x ) 

Assign a 1; ..., a t to interval coefficients and name it S'*; 

-< a := an arbitrary monomial ordering on a\,.. . , a t , 

E,N-.= { },{!}; 

-< x ,a : =The block ordering of -< x , -< a 

Return PGB -MAIN(P, E, A, -< x a , L); \\L is the ordered set of inter¬ 
val coefficients which is needed to check consistency 

end procedure 


The PGB-MAIN algorithm is the same which which is used in PGB 
algorithm. We only change the definition of consistency as below. 

Definition 6.12. Let [an, Pi),..., [a*, (3 t ) be t real intervals and E,N C 
R[ai,..., a t ]. The pair ( E, N ) is called consistent if it is not redundant, or 
equivalently, 


[V(E) \ V (N)] n [«!, A) x • • • x [at, Pt) ^ 0. 

According to the above definition, we change the IsConsistent algo¬ 
rithm to Interval-IsConsistent, which checks the consistency for radical 
membership and redundancy determination. 


23 





Algorithm 5 INTERVAL-IsCONSISTENT 

procedure lNTERVAL-IsCoNSlSTENT(i?, N, [« 1 , /3i), ..., [a t , fit)) 
test :=false; 
flag :=false; 

if (E, N) is not redundant then 
test :=true; 

end if 
if test then 
flag :=false; 

for g E N while flag = false do 
if g f \J (E) then 
flag :=true; 

end if 
end for 
end if 

Return flag ; 
end procedure 

Remark 6.13. It is worth noting that redundant triples will be omitted before 
the algorithm goes to continue with them. This property causes that IGS 
returns less triples than PGB. 

7 Examples 

7.1 Solving Fuzzy Polynomial Systems 

In this section we state the ability of interval Grobner system to solve para¬ 
metric polynomial systems for which the parameters range over specific in¬ 
tervals, which may appear, for instance, when analysing fuzzy polynomial 
systems. In doing so, we can use the same IGS algorithm, after converting 
an interval polynomial system to a parametric one. To give an example, we 
point at resolution of fuzzy polynomial systems which converts to solve a 
parametric polynomial system with parameters range over [0,1] (see [1]). 

Example 7.1. ([1, Example f.3] resolved) Consider the following system of 
fuzzy polynomials: 

f (0.25,0.75,0.625) = (0, 0.5, 0, 5)x 5 + (0.5, 0.5, 0.25)y 
{ (1.25,1.25,0.75) = (0.5,1, 0.5)a: 4 + (1.5, 0.5, 0.5)xy 
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To solve this system, the general way is to decompose it into four paramet¬ 
ric polynomial systems. Here we solve only the third and fourth systems. 
Consider the third system as follows: 

[ h = (| - \h)x b + {\h)y + | - \h 

/o\ . J f 2 = + \h)x b + (| - \h)y - | + §/i 

‘ | h = (-§ + h)x 4 + (2 - \h)xy - f h 

( U = (! - \h)x A + (1 + \h)xy - 2 + \h 

Computing an interval Grobner basis w.r.t. y -<i ex x, we have only a triple 
({ }; {h ~ 1}, {1}) which means that we have no solutions whenever h ^ 
1. As our algorithm deals with the values of h e [0,1), we must consider 
the case h — 1 seperately. hi this case, we find a solution (h — l,x — 
— 1.473157368,?/ = 0.5) regarding the sign of variables. 

The fourth system is: 

[ h = (| - \h)x b + (f - \h)y + | - \h 

(a) . J h = (-5 + \ h ) xb + (\ h )y -1 + \ h 

I h = (~| + /i)x 4 + (1 + - |/i 

[ ,/ 4 = (1 - + (2 - i/i)xj/ - 2 + |/i 

where h e [0,1], x < 0 and y < 0. Computing an interval Grobner basis, we 
see only one triple ({ }, {/? 2 — 2/z + 3}, {1}) for which its Grobner basis equals 
to {1}. Therefore the above system has no solution for h G [0,1). Note that 
our method deals with [0,1) here and so we must check the system for h — 1 
separately. By this, the system has no solution again by considering the sign 
of x and y. and so this system has no solution. 

Note that another way to inform that these system have or have not any 
real solutions (and not the exact form of solutions), is to check their con¬ 
sistency by what we have stated in Algorithm 6.1, since all of variables have 
interval form in these systems. In doing so we must add three auxiliary poly¬ 
nomials h+(h—l)h 2 , x-\-x 2 and y—y 2 to the System (3) and h+(h—l)h 2 , x+x 2 
and y + y 2 to the System (4), where h,x and y are some extra real variables. 
The result of course is that (3) is consistent however (4) is not. 

7.2 Real Factors 

One of the interesting problems in the context of interval polynomials is the 
Divisibility Problem stated as follows (See [32]): 


25 



For an interval polynomial [/] G [K]{xi,... ,x n } and a real polynomial 
g G M[xi...., x n ], determine whether there is a polynomial p G .F([/]) such 
that g is a factor of p. 

In [32] there is stated an efficient method based on linear programming 
techniques and nice properties of polytopes. In the sequel we exalin our 
method to solve this problem by interval Grobner basis. To continue, let 
us say that g i-divides [/] whenever the answer of the above problem is yes 
(Note that the letter ”i” stands for interval). By this conception we can now 
explain the following criterion based on interval Grobner basis. 

Theorem 7.2. Using the above notations, let [S] = {[/], g} C [M]{aq,..., x n } 
and [Q] be a reduced interval Grobner system for [S'] with respect to a mono¬ 
mial ordering Then g i-divides [/] if and only if there exists a triple 
(E, N, G ) G [Q] where G = {g} and V(E) \ V(N) ± {0}. 

Proof. It is easy to see that g i-divides [/] if and only if there exists an 
specialization a for which g\a([f]) and of course <r([/]) fz 0, by the statement 
of divisibility problem. This implies that uQS 1 ]) can be expressed only by 
{q | and therefore there exists a pair of parametric sets (E.N) such that 
(E,N,{g})t={g]. □ 

Example 7.3. We are going to solve the divisibility problem for 

[/] = [-1, l}x 2 + [-3,1 }y 2 + [1/2, 2 }xy G [M]{x, y} 

and g = x — \[2y G M[x,t/]. By computing a reduced interval Grobner basis 
for [S] = {[f],g}, we find the triple 

({cV2 + 2a + b},{l},{x-V2y}) 

where a,b and c denote inner values of the intervals [—1,1], [—3,1] and 
[1/2,2], This means that if the values of a,b and c satisfy the equation 
c \/2 + 2a + b = 0, then there exists some p G E'df]) such that g\p. For in¬ 
stance, by evaluating a — 1/3, b = —2 and so c = 2\/2/3 we find p = l/3x 2 — 
2y 2 +2y/2/3xy G lE([f]) which is divided by g (Note thatp = l/3(x+3y/2y)g). 
Therefore g i-divides [/]. 

Remark 7.4. We can use Theorem 7.2 for further aims. Let f, g G M[aq,..., x n ] 
where g does not divide f. Then one can use Theorem 7.2 to find a polyno¬ 
mial f with the same coefficients of f which contain a few perturbation and 
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g | f. In doing so, one can convert f to an interval polynomial [/] by putting 
the interval [c — e, c + e] instead of the coefficient c, for each coefficient c 
appearing in f. Then using Theorem 7.2 one can increase the e up enough 
until g i-divides [/] with the desired precision. 

8 Conclusion 

In the current paper we have introduced the concept of interval Grobner sys¬ 
tem as a novel computational tool to analyse interval polynomial systems. 
We have further designed a complete algorithm to compute it using the ex¬ 
isting methods to analyse parametric polynomial systems. This concept can 
solve some important problems from the applied and engineering research 
fields such as solving fuzzy polynomial system, as devoted in this paper. 
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